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Abstract 

We predict the branching ratios of Z — > e^/i 1 * 1 , Z — > e^T^ and Z — ► y^r decays in 
the model III version of the two Higgs doublet model, with the inclusion of one and two 
spatial non-universal extra dimensions. We observe that the the branching ratios are not 
sensitive to a single extra dimension, however, this sensitivity is considerably large for two 
extra dimensions 
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1 Introduction 

The lepton flavor violating (LFV) interactions are interesting in the sense that they are sensitive 
the physics beyond the standard model (SM) and they ensure considerable information about 
the restrictions of the free parameters, appearing in the new models, with the help of the 
possible accurate measurements. Among LFV interactions, the Z decays with different lepton 
flavor outputs, such as Z — > e/i, Z — > er and Z — > fir, are rich enough to study and there is an 
extensive work related to these decays in the literature PQ-[I2]- The Giga-Z option of the Tesla 
project which aims to increase the production of Z bosons at resonance ^3] stimulates one to 
make theoretical works on such Z decays. 

In the framework of the SM the lepton flavor is conserved and, for the flavor violation in the 
lepton sector, there is a need to extend the SM. One of the candidate model is so called z/SM, 
which is constructed by taking neutrinos massive and permitting the lepton mixing mechanism 
[T4*] . In this model, the theoretical predictions for the branching ratios (BRs) of the LFV Z 
decays are extremely small in the case of internal light neutrinos PjJ |2] 

BR(Z -> e ± /i ± ) ~ BR(Z -> e ± r ± ) ~ 10" 54 , 

BR(Z -> /i ± r ± ) < 4 x KT 60 . (1) 

They are far from the experimental limits obtained at LEP 1 [TK] : 
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and from the improved ones at Giga-Z 6\: 
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with / = 0.2 — 1.0. Notice that these numbers are obtained for the decays Z — ► IJ2 + hh, 
namely 

BR(Z -> If If) = nz^hk + kh) (4) 

t z 

To enhance the BRs of the corresponding LFV Z decays some other scenarios have been studied. 
The possible scenarios are the extension of z/SM with one heavy ordinary Dirac neutrino 
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the extension of z/SM with two heavy right-handed singlet Majorana neutrinos [2], the Zee 
model [Zj, the model III version of the two Higgs doublet model (2HDM), which is the minimal 
extension of the SM jH], the supersymmetric models [HUH]], top-color assisted technicolor model 

cu. 

The present work is devoted to predictions of the BRs of Z — > e // , Z — > e ± T ± and 
Z — > /i ± r ± decays in the model III version of the 2HDM, with the inclusion of one and two 
spatial extra dimensions. Our motivation is to check whether there is an enhancement in the 
BRs of these decays due to the extra dimensions. The possible existence of new dimensions 
reach great interest recently and there is a large amount of work done in the literature [16 - [32J . 
The idea of extra dimensions was originated from the study of Kaluza-Klein j33] which was 
related to the unification of electromagnetism and the gravity and the motivation increased 
with the study on the string theory which was formulated in a space-time of more than four 
dimensions. Since the extra dimensions are hidden to the experiments at present (for example 
see [30 J, the most favorable description is that these new dimensions are compactified to the 
surfaces with small radii, which is a typical size of corresponding extra dimension. This leads 
to appear new particles, namely Kaluza-Klein (KK) modes of the particles in the theory. In the 
case that all the fields feel the extra dimensions, so called universal extra dimensions (UED), the 
extra dimensional momentum, and therefore, the KK number at each vertex is conserved. The 
compactification size R has been predicted as large as few hundereds of GeV [T71 [TBI EH l2*U] . 
in the range 200 — 500 GeV, using electroweak precision measurements [21], the B — B -mixing 
[22[ 123] and the flavor changing process b —>■ sj [21] for a single UED. 

The assumption that the extra dimensions are at the order of submilimeter distance, for two 
extra dimensions, the hierarchy problem in the fundamental scales could be solved and the true 
scale of quantum gravity would be no more the Planck scale but in the order of electroweak 
(EW) scale ^21 El- in this case, the gravity is spreading over all the volume including the extra 
dimensions, however, the matter fields are restricted in four dimensions, called four dimensional 
(4D) brane, or in 4D surface which has a non-zero thickness in the new dimensions, called fat 
brane (see for example [25j). This type of extra dimensions, accessible to some fields but not all 
in the theory, are called non-universal extra dimensions (NUED). Contrary to the UED, in the 
NUED, the KK number at each vertex is not conserved and tree level interaction of KK modes 
with the ordinary particles can exist. The study in j2H] is devoted to the appearence of a very 
light left handed neutrino in the NUED where only the right handed neutrino is accessible to 
the extra dimension. In the another work [2Zj, the effect of brane kinetic terms for bulk scalars, 
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fermions and gauge fields in higher dimensional theories, have been studied. In [28] the electric 
dipole moments of fermions and some LFV decays have been analyzed in the framework of 
NUED. 

Here, we predict the BRs of the LFV Z decays in the model III with the assumption that 
the extra dimensions are felt by the new Higgs doublet and the gauge sector of the theory. The 
Z decays under consideration are induced by the internal neutral Higgs bosons h° and A and 
their KK modes carry all the information about the new dimensions, after the compactification 
of the single (double) extra dimension on the orbifold S 1 / Z2 ( (5* 1 x S 1 )/Z2- Here, the KK 
number in the vertices is not conserved, in contrast to the UED case. The non-zero KK modes 
of neutral Higgs fields H have masses J m 2 H + m 2 n (\Jm 2 H + m 2 n + m 2 ) with m& = k/R in one 
(two) extra dimension. We observe that the BRs of the processes we study enhance almost two 
orders larger compared to the ones without the extra dimensions, in the case of two NUED, 
since the neutral Higgs KK modes are considerably crowded. 

The paper is organized as follows: In Section 2, we present the effective vertex and the BRs 
of LFV Z decays in the model III version of the 2HDM with the inclusion of NUED. Section 3 is 
devoted to discussion and our conclusions. In appendix section, we give the explicit expressions 
of the form factors appearing in the effective vertex. 

2 Z —> l^l} decay in the model III with the inclusion of 
non-universal extra dimensions. 

The extension of the Higgs sector in the SM brings new contributions to the BRs of the processes 
and makes it possible to obtain the flavor changing neutral current (FCNC) at tree level, which 
plays an important role in the existence of flavor violating (FV) interactions. Therefore, the 
multi Higgs doublet models are worthwhile to study. The 2HDM is one of the candidate for 
the multi Higgs doublet models. In the model I and II versions of the 2HDM, the FCNC at 
tree level is forbidden, however, those type of interactions are possible in the model III version 
of the 2HDM. The lepton flavor violating (LFV) Z decay Z — > l^l^ can be induced at least in 
the one loop level in the framework of the model III. 

The addition of possible NUED, which are experienced by the gauge bosons and the new 
Higgs particles, brings new contributions to the BRs of the decays under consideration. In 
the model III, the part of lagrangian which carries the interaction, responsible for the LFV 
processes in 5 (6) dimension, reads 

£y = f§ 6 )i/i (&li/(*)=o) Ej + h.c. , (5) 
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where the couplings £5(6) ij are 5(6)-dimensional dimensionful Yukawa couplings which induce 
the LFV interactions. These couplings can be rescaled to the ones in 4-dimension as £5(6) ij = 
\J2-kR (2irR) with lepton family indices i,j 1 . Here, 2 is the new scalar doublet, R is the 
compactification radius, k and Ej are lepton doublets and singlets, respectively. The scalar and 
lepton doublets are the functions of and y (y, z), where y (z) is the coordinate represents 
the 5(6) 'th dimension. Here we assume that the Higgs doublet lying in the 4 dimensional brane 
has a non-zero vacuum expectation value to ensure the ordinary masses of the gauge fields and 
the fermions, however, the second doublet, which is accessible to the extra dimensions, has no 
vacuum expectation value, namely, we choose the doublets 0i and 2 and the their vacuum 
expectation values as 



\ , ( V2x 
v + H° ) + 1 i X ° 



and 



< <h >= ^ ( I ) ; < 02 >= . (7) 

This choice ensures that the mixing between neutral scalar Higgs bosons is switched off and it 
would be possible to separate the particle spectrum so that the SM particles are collected in 
the first doublet and the new particles in the second one 2 . Here we consider the gauge and 
CP invariant Higgs potential in two extra dimensions 

\/(0!,0 2 ) = 6(y)6(z)c 1 (<i>t<i>i-v 2 /2) 2 + c 2 (<j>t<h) 2 + 

+ 5{y) 8{z) (c 3 [(0M - v 2 /2) (0+0 2 )] + c^tfaMfo) 

- (0+0 2 )(0+0O] + c 5 [ J Re(0+0 2 )] 2 + c 6 [/m(0+0 2 )] 2 ) + c 7 , (8) 

with constants q, i = 1, 7. 

Since, only the new Higgs field 2 is accessible to extra dimensions in the Higgs sector, there 
appear KK modes 2 n ' r ' ) of 2 in two spatial extra dimensions after the compactification of the 
external dimensions on an orbifold (S 1 x S 1 )/Z 2 , 

1 f 00 1 

M*,V,z) = ^^|0r ) (^) + 2 d/2 E# r) (^) co ^/^ + ^)) , (9) 

where d = 2, the indices n and r are positive integers including zero, but both are not zero at 
the same time. Here, <p2'°\ x ) is the 4-dimensional Higgs doublet which includes the charged 

1 In the following, we replace £ D with £jy where "N" denotes the word "neutral" 
2 Here H 1 (H 2 ) is the well known mass eigenstate h a (A ). 



4 



Higgs boson H + , the neutral CP even (odd) Higgs bosons h° (A ). The KK modes of the 
charged Higgs boson (neutral CP even (odd) Higgs h° (A )) have the mass J ni 2 H± + m\ + m 2 
( J m\ + + m 2 {yJm 2 A o + ™ 2 n + rn%) ), where m n = n/R and m r = r/R. Furthermore, we 
assume that the compactification radius R is the same for both new dimensions. Notice that 
the expansion for a single extra dimension can be obtained by setting d — 1, taking z — 0, and 
dropping the summation over r. In addition to the new Higgs doublet, also the gauge fields feel 
the extra dimensions, however, there is no additional contribution coming from the KK modes 
of Z boson in the process under consideration since the Z boson does not enter into calculations 
as an internal line. The Z — h° KK mode-A° KK mode vertex is the same as the 4-dimensional 
one after the integration over extra dimensions, except a small correction of the coupling due 
to the gauge field 0-mode-KK mode mixing (see section 6 for details). 

Now, we would like to present the general effective vertex for the interaction of on-shell 
Z-boson with a fermionic current: 

r„ = 7 „(/v - fA 75) + — {fu + f E 75) <v„ q v (10) 

where q is the momentum transfer, q 2 = (p — p') 2 , fy (Ja) is vector (axial-vector) coupling, 
/m (/e) magnetic (electric) transitions of unlike fermions. Here p (— p') is the four momentum 
vector of lepton (anti-lepton) (see Fig. El for the necessary 1-loop diagrams due to neutral 
Higgs particles). Since the LFV Z boson decay exists at least in the loop level, the KK modes 
of neutral Higgs particles h° and A contribute to the self energy and vertex diagrams as internal 
lines. The leptons live in the 4D brane and therefore they do not have any KK modes. Notice 
that, in the case of non- universal extra dimension, the KK number needs not to be conserved 
and there exist lepton — lepton — h° KK mode (A KK mode) vertices which can involve two 
zero modes and one KK mode. 

The vector (axial-vector) fy (/a) couplings and the magnetic (electric) transitions fu (/b) 
including the contributions coming from a single extra dimension can be obtained as 

00 

fy = /«?> + 2 5;/<"\ 

n=l 



>3C 



f A = /i 0) + 2^#, 

n=l 

00 

JM - J M + 1 2^ JM > 

n=l 

00 

fE = ff+^fP, (11) 
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where fy A M E are the couplings in the 4-dimensions and fy\ M E are the ones due to the KK 

fin) 

1 V,A,M,. 

,(0) 



modes of the scalar bosons S = h°,A°. The KK mode contributions fy~AME can be easily 



obtained by replacing the mass squares m s in fvAME by Trig + m^, with m n = n/R and the 
compactification radius R. We present the explicit expressions for the couplings fy\ M E in 
the appendix, by taking into account all the masses of internal leptons and external lepton 
(anti-lepton). 

If we consider two NUED, the couplings fyA me appearing in eq. should be replaced 
by fyA me an d the summation would be done over n,r = 0, 1,2... except n = r = 0. Here 
fvAM e can be obtained by replacing the mass squares m| in fy A M E by m| + m\ + m^, with 
TTin — n/R, m r = r/R. Furthermore, the number 2 in front of the summations in eq. (jllj) 
would be replaced by 4. 

Finally, the BR for Z — > IT If can be written in terms of the couplings fy, fA, f\i and 

as 

BH(Z ir i2+ ) = J- ^ fifyf + ]!a? + ^(I/mI 2 + l/ £ | 2 )} (12) 

2 

where aw = f- and Y z is the total decay width of Z boson. In our numerical analysis we 
consider the BR due to the production of sum of charged states, namely 

BR(Z -> If If) = T (Z^(kk + l2li) (13) 

3 Discussion 

The LFV Z decays Z — > e fi , Z — > e ± r ± and Z — > /i ± r ± strongly depend on the Yukawa 
couplings ^jj 3 ,i,j = e, /x, r in the model III version of 2HDM and these couplings are free 
parameters which should be restricted by using the present and forthcoming experiments. At 
first, we assume that the couplings which contain at least one r index are dominant similar to 
the Cheng-Sher scenario [33] and, therefore, we consider only the internal r lepton case among 

D 

N,ij 

N,rei S7V,T/i 

N,tt- 

The upper limit of Cat,™ i s predicted as 30 GeV (see [33] and references therein) by using the 
experimental uncertainty, 10~ 9 , in the measurement of the muon anomalous magnetic moment 
and assuming that the new physics effects can not exceed this uncertainty. Using this upper 



others. Furthermore, we assume that the Yukawa couplings are symmetric with respect 
to the indices i and j. As a result, we need the numerical values for the couplings £^ re , £ D 
and fjP 



3 The dimensionfull Yukawa couplings are defined as = y ^f" ^N,ij 



limit and the experimental upper bound of BR of fi — > ey decay, BR < 1.2 x 10~ n , the coupling 
£® Te can be restricted in the range, 1CT 3 — 1CT 2 GeV (see [US]). For the Yukawa coupling £® TT , 
we have no explicit restriction region and we use the numerical values which are greater than 
£jvr/r Furthermore, the addition of the extra dimensions bring new parameter, namely the 
compactification radius R which arises from the compactification of the a single (double) extra 
dimension on the orbifold S 1 /Z 2 ((S 1 x S 1 )/Z 2 ). 

In the present work, we study the prediction of the NUED on the BR of the LFV processes 
Z — > if if, in the framework of the type III 2HDM. We see that the contribution coming from 
two extra dimensions are considerably large compared to the one coming from a single extra 
dimension, due to the crowd of neutral scalar Higgs boson KK modes. 

Throughout our calculations we use the input values given in Table ((H). 



Parameter 


Value 


m T 

mw 
m z 
Gp 
T z 

sin 9\y 


0.106 (GeV) 

1.78 (GeV) 

80.26 (GeV) 

91.19 (GeV) 

1.1663710- 5 (Ge1/- 2 ) 

2.490 (GeV) 

V0.2325 



Table 1: The values of the input parameters used in the numerical calculations. 

Fig. El is devoted to £j} jTe dependence of the BR (Z — > ^ e*) for ^ jT „ = 1 GeV, m h o = 
100 GeV and m^o = 200 GeV. The solid-dashed-small dashed lines represent the BR without 
extra dimension-including a single extra dimension for 1 / R = 500 Ge V-including two extra 
dimensions for 1/R = 500 GeV. It is observed that the BR is not sensitive to the extra 
dimension effects for a single extra dimension. However, for two NUED, there is a considerable 
enhancement, almost two orders, in the BR compared to the one without extra dimensions, even 
for the small values of the coupling £jy T „. This is due to the crowd of neutral Higgs boson KK 
modes. This enhancement can be observed also in Fig. 0] where the BR is plotted with respect 
to the compactification scale 1/R for ^Nre = 0.05 GeV, £j} T „ = I GeV, m h o = 100 GeV and 
m^o = 200 GeV. In this figure the solid-dashed-small dashed lines represent the BR without 
extra dimension-including a single extra dimension- including two extra dimensions. It is seen 
that in the case of two extra dimensions the BR reaches almost twice of the one without extra 
dimensions, for the values of the compactification scale, 1/R — 2000(^6^. This enhancement 
becomes negligible for the larger values of the compactification scales, 1/R > 5000 GeV. The 
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possible enhancement due to the effect of two NUED on the theoretical value of the BR of the 
corresponding Z decay is worthwhile to study. 

In Fig. we present £j} TT dependence of the BR(Z — > e^) for £jy Te = 0.05 GeV, 
m^o = 100 GeV and m^o = 200 GeV. The solid-dashed-small dashed lines represent the BR 
without extra dimension-including a single extra dimension for 1/R — 500 GeVMncluding two 
extra dimensions for 1/R — 500 GeV. Similar to the previous process, the BR is not sensitive 
to the extra dimension effects for a single extra dimension and this sensitivity increases for two 
NUED. The enhancement of the BR of two NUED is more than two orders larger compared 
to the one without extra dimensions. Fig. |H1 is devoted to the compactification scale 1/R 
dependence of the BR for £§ >TT = 10 GeV, £§ >Tfl = 1 GeV, m h o = 100 GeV and m A o = 200 GeV. 
In this figure the solid-dashed-small dashed lines represent the BR without extra dimension- 
including a single extra dimension- including two extra dimensions. The enhancement in the 
BR for the intermediate values of the compactification scale, namely 1/R ~ 1000 GeV, is more 
than one order. Similar to the previous decay, this enhancement becomes small for the larger 
values of the compactification scales, 1/R> 5000 GeV. 

Finally, Fig. (JBJ) is devoted to the <^ TT (the compactification scale 1/R) dependence 
of the BR of the decay Z -»• r*/^ for ^, r/1 = 1 GeV (^ >TfM = 1 GeV, £% jTT = 10 GeV), 
m h o = 100 GeV and m^o = 200 GeV. In Fig. [7|the solid-dashed-small dashed lines represent the 
BR without extra dimension- including a single extra dimension for 1/R = 500 Gel^-including 
two extra dimensions for 1/R = 500 GeV. In the case of two extra dimensions, even for 
small Yukawa couplings, it is possible to reach the experimental upper limit of the BR of the 
corresponding decay, since the enhancement in the BR is two order larger compared to the 
case without extra dimensions. In Fig. |Hl the solid-dashed-small dashed lines represent the BR 
without extra dimension-including a single extra dimension-including two extra dimensions. It 
is observed that, in the case of two extra dimensions, the BR reaches almost twice of the one 
without extra dimensions, even for intermediate values of the compactification scale, 1/R = 
2000 GeV. For the larger values of the compactification scales, 1/R > 5000(^6^, there is no 
enhancement in the BR of the present decay. 

At this stage we would like to present our results briefly. 

• With the inclusion of a single NUED, the enhancement in the BR of the LFV Z decays 
is small for the intermediate values of the compactification scale 1/R. 

• In the case of two NUED, even for the small values of the Yukawa couplings, it is possible 
to reach the experimental upper limits of the BRs of the LFV Z decays, since the en- 
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hancement in the BR is two order larger compared to the case without extra dimensions 
for the intermediate values of the compactification scale 1/R. This enhancement is due to 
the crowd of the KK modes and it is an interesting result which may ensure an important 
information to test the existence of the NUED, and if it exists, to decide its number and 
to predict the lower limit of the compactification scale, with the help of more accurate 
experimental results. 

As a summary, the effect of two NUED on the BRs of LFV Z decays Z—tlfl^ is strong and 
the more accurate future experimental results of these decays will be useful to test the possible 
signals coming from the extra dimensions. 
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5 The explicit expressions appearing in the text 



Here we present the explicit expressions for f v , and f% |H] (see eq. (fTTjO : 
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The parameters cy and c A are = — | and cy = 4 — sin 2 0^. In eq. (fTHjl the flavor changing 
couplings C^iji represent the effective interaction between the internal lepton i, (i = e, r) 
and outgoing (incoming) j = 1 (j = 2) one. Here the couplings ^,1^ are complex in general 
and they can be parametrized as 



D I i^- 
NAL I e 



(17) 



where i, lj denote the lepton flavors and % are CP violating parameters which are the possible 
sources of the lepton EDM. However, in the present work we take these couplings real. 
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6 Gauge boson mass matrix and gauge coupling 



Here we study an abelian model in the case of a single extra dimension (two extra dimensions) 
with two Higgs fields, where one of Higgs field, 4>i(x), is localized at the y = (y = z = 
0) boundary of the S l /Z 2 ((S 1 x S 1 )/Z 2 ) orbifold and the other one (f>i(x,y) (4>i(x, y, z)) is 
accessible to the extra dimension(dimensions). Furthermore, we choose that only the first 
Higgs field has a non-zero vacuum expectation value, and, including a single extra dimension, 
the Higgs fields read: 

4> 2 (x) = -^=(h 2 (x,y) +ix2(x,yj) • (18) 

Our aim is to obtain the gauge boson mass matrix, which is not diagonal in the case of the 
non-universal extra dimensions where the the gauge sector and the Higgs field <fi 2 is accessible to 
the extra dimensions but the first Higgs field does not. In the case of a single extra dimension, 
the detailed analysis on this issue has been done in [3T] and both Higgs fields assumed to have 
vacuum expectation values in that work. We will present the crucial steps of this work briefly 
and we repeat the same analysis for two extra dimensions. 

The part of Lagrangian which carries the gauge and Higgs sector in a single extra dimension 
reads 

£(x,y) = -- A F MN F MN + {D M< p 2 Y{D M <p 2 )+5{y){D ll< p i r{D^ l ) 

-V(fa,<h) + C GF (x,y), (19) 

where V is the CP and gauge invariant Higgs potential, D M = 8m + ie 5 A M (x, y), (M = /i, 5) 
is the covariant derivative in 5 dimension and CcF{x,y) is the gauge fixing term. 

Now, we will present the sources of the gauge boson mass matrix in the lagrangian eq. (jl9j) : 

1. F b »F^ in the part F MN F MN , 

2. (A^A") in the part (£> M 0i)* (£>"0i) , 

where F MN = d A N — d M A , and, after compactification, the gauge fields A N read 

M*>v) = ^^{4°H^) + 2 1/2 E4 n) ^)cos(^/ J R)| , 

As(x,y) = ^y^|2 1/2 E4 n) (^)sin(^/ J R)| . (20) 
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Notice that the (-D M 2 )* (D M <p 2 ) term in the lagrangian (see eq. ()19j)) does not produce any- 
mass term for the gauge field since the scalar field 02 does not have any vacuum expectation 
value. The integration over the extra dimension y results in the mixing of zero mode and KK 
mode gauge fields and the gauge boson mass matrix is obtained as (see [HI]): 



M 2 



\/2m 2 



V2m 2 



( m 2 

^2m 2 2 m 2 + (l/R) 2 2 m 2 

1 2 m 2 2m 2 + (2/R) 2 



\ 



^2m~ 



(2l\ 



1 



V : : : 

where m 2 = e 2 v 2 . Using the determinant equation 

oo oo 

det (M 2 - A /) = ( J] (n 2 /R 2 - A)) (m 2 - A - 2 A m 2 £ ~ {n/ R)2 _ A 
the eigenvalues of the matrix is obtained by solving the transcendental equation 

m A („) = 7r m 2 R cot(7r m A ( n ) R) , 



(22) 



(23) 



and KK mass eigenstates are given by 



\ -1/2 oo 

l + vr 2 m 2 J R 2 + ^ 



- 2m AW m / 1 

L m 2 _ ij/R) 2{ V z) A M) 



j=0 m A(") 



(24) 



For the non-abelian case the gauge field mass spectrum is analogous to the abelian one presented 
above and the transcendental equation for Z boson is 



m Z ( n ) = -Km 2 z R cot(n m Z (n) R) 



(25) 



and the corresponding coupling reads 



few 



V2g 1 + 



m 



+ 



7T 



2 R 2 m%\- 1 / 2 



rrir 



At this stage we try to make the same analysis for two extra dimensions. The part of Lagrangian 
which carries the gauge and Higgs sector in two extra dimensions is 

1 



(26) 



C(x,y) = -^F MN F MN + (D M( j )2 r(D M ( j )2 ) + 5(y)5(z)(D ll( j )i r(D^ l ) 
V(<j)i, 2 ) + C GF (x,y,z) , 



(27) 



Dm = 9m + ie&AM(x,y, z), (M = //, 5, 6) is the covariant derivative in 6 dimension. In this 
case the sources of the gauge boson mass matrix in the lagrangian (eq. ([27))) are: 
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1. F 5 ^ F 5fl and F 6/i F 6/i in the part F MN F MN 



2. (A^A^) in the part (D^cfti)* (D^cfri), and, after compactification, the gauge fields A N read 

A,(x,y) = j^L^°\x)+2f:A^\x)co8(ny/R + rz/R)\ , 

1 f 00 1 

A 5i 6)(x,y) = j^hJ2A^(x)sin(ny/R + rz/R)\ , (28) 

The integration over the extra dimensions y and z results in the mixing of zero mode and KK 
mode gauge fields similar to the one extra dimension case and the gauge boson mass matrix is 
obtained as: 




(29) 

with m 2 = e 2 v 2 . Here, the mass spectrum is richer compared to a single extra dimension case. 
Now the determinant equation reads 

det (Mi - A J) = - A) ) (m 2 - A - 4Am 2 E r ) = , (30) 

where the indices n and r are all positive integers including zero, but both is not zero at the 

same time, and the transcendental equation to obtain the eigenvalues of the matrix is 

2 m 2 1 

m%) = -(-^r) 2 + + 10 m 2 ) m 2 A(p) - (m A(p) m 2 n R') (— - m 2 A(p) ) cot (7rm A(p) R') 

- (m 2 A(p) m 2 ) (— -m\ {p) ) ^ ^= cot \/K R') , (31) 

where A r = m 2 A(p) — r 2 /R' 2 , R' = R/V2 and p is the positive integer. Similar to the one 
extra dimension, for the non-abelian case, the gauge field mass spectrum is analogous to the 
abelian one presented above and the transcendental equation is obtained by replacing the mass 
m A („) in eq. (}3*Tj) by m Z (n). For two extra dimensions there appears a new gauge coupling 
due to the complicated mass mixing, however, in our numerical calculations, we used the one 
obtained in the single extra dimension case by expecting that the new contributions do not 
affect the behaviors of the physical parameters we study. Notice that this coupling enters in 
the calculations only for the zero mode Z boson case, since there is no diagram which includes 
the KK mode virtual Z bosons. 
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7 The vertices appearing in the present work 



In this section we present the vertices appearing in our calculations. For a single extra dimension 
the Z boson gauge coupling is given in eq. (j2Bj) and, in the present work, this coupling reads 



where m Z (o) is obtained by solving the eq. (|25j) for n = 0. Furthermore, since there is no mixing 
between neutral scalar Higgs bosons H° and h° due to the our choice (see the section 2), the 
tree level interaction — H° — A does not exist. Here L and R denote chiral projections 
L(R) = 1/2(1 =)= 75), the parameters cl(r) read, cl = — 1/2 + s^, cr = and cw = cos 0w, 
sw = sin 9w, where 6w is the weak angle. 




(32) 
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(n, r) 



(n, r) 



-igz 



^ [c L L + c R R] 



9 ±M + viY 



lew 



2V2 



>J1 



1 

2V2 



- ) + m + & )75 



D 



5 y 



Figure 1: The vertices used in the present work. 
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(a) 




(c) (d) 



Figure 2: One loop diagrams contribute to Z — > k + j~ decay due to the neutral Higgs bosons 
ho and A in the 2HDM. % represents the internal, j (k) outgoing (incoming) lepton, dashed 
lines the vector field Z, h and A fields. In the case 5 (6) dimensions the vertices are the same 
but there are additional contributions due to the KK modes of h and A fields. 
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Figure 3: f£ jTe dependence of the BR(Z -> /i ± e ± ) for £% >Tli = 1 GeV, m h o = 100 GeV 
and rriA° = 200 GeV. The solid-dashed-small dashed lines represent the BR without extra 
dimension-including a single extra dimension for 1/R = 500 GeV^-including two extra dimen- 
sions for l/R = 500 GeV. 
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Figure 4: The compactification scale 1/R dependence of BR (Z — > ^ e ± ) for £jv re = 0.05 GeV, 
^N,Tfi = ^GeV, m h o = 100 GeV and m A o = 200 GeV. The solid-dashed-small dashed lines 
represent the BR without extra dimension-including a single extra dimension-including two 
extra dimensions. 
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Figure 5: £,n,tt dependence of the BR (Z -> r ± e ± ) for ^ Te = 0.05 GeV, m h o = 100 GeV 
and mjsfi = 200 GeV. The solid-dashed-small dashed lines represent the BR without extra 
dimension-including a single extra dimension for 1/R = 500 GeV-including two extra dimen- 
sions for 1/R = 500 GeV. 
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Figure 6: The compactification scale 1/R dependence of the BR(Z — > r ± e ± ) for £® TT = 
lOGeV, £% jTli = 1 GeV, m h o = 100 GeV and m A o = 200 GeV. The solid-dashed-small dashed 
lines represent the BR without extra dimension-including a single extra dimension-including 
two extra dimensions. 



25 



i I I I I I I I I I I 

'5 10 15 20 25 30 35 40 45 50 



Figure 7: The £§ TT dependence of the BR of the two decay Z — > fj^ for £§ Tfi = 1 GeV, 
rriho = 100 GeV and m^o = 200 GeV. The solid-dashed-small dashed lines represent the BR 
without extra dimension-including a single extra dimension for 1/R — 500 GeV-including two 
extra dimensions for 1/R — 500 GeV. 
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Figure 8: The compactification scale 1/R dependence of the BR of the two decay Z — > pi for 
^N,rfj, = I GeV, ^ )TT = WGeV, m h o = 100 GeV and m A o = 200 GeV. The solid-dashed-small 
dashed lines represent the BR without extra dimension-including a single extra dimension- 
including two extra dimensions. 
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